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ABSTRACT 



Context. The connection of cluster mergers with the presence of extended, diffuse radio sources in galaxy clusters is still being 
debated. 

Aims. We aim to obtain new insights into the internal dynamics of Abell 959, showing evidence of a diffuse radio source, analyzing 
(-H ' velocities and positions of member galaxies. 

Methods. Our analysis is based on redshift data for 107 galaxies in the cluster field acquired at the Telescopio Nazionale Galileo. We 
also use photometric data from the Sloan Digital Sky Survey (Data Release 6). We combine galaxy velocities and positions to select 81 
galaxies recognized as cluster members and determine global dynamical properties. We analyze the cluster searching for substructures 
by using the weighted gap analysis, the KMM method and the Dressler-Shectman statistics. We also study the 2D galaxy distribution 
CZ5 , in the field of the cluster. We compare our results with those from X-ray and gravitational lensing analyses. 

Results. We estimate a cluster redshift of (z) = 0.2883 ± 0.0004. We detect an NE high velocity group at 5' from the cluster center 
with a relative line-of-sight (LOS) velocity of ~ +1900 km s _1 with respect to the main system. We also detect a central, dense 
structure elongated along the SE-NW direction likely connected with the two dominant galaxies and their surrounding cores. This 
elongated central structure is probably the trace of an old cluster merger. The LOS velocity dispersion of galaxies is very high. By 



excluding the NE clump we obtain tr v = 1025*^ km s" 1 and a virial mass M(< R = 1.48 ft^jMpc) = 1.15^ xlO 15 M Q . Our 
results suggest that this cluster is forming along two main directions of mass accretion. 

Conclusions. Abell 959 is confirmed to show the typical characteristics of radio clusters; i.e., it is very massive and shows a young 
dynamical state. However, deeper radio observations are needed to clarify the nature of the diffuse radio emission in Abell 959 and 
its connection with the cluster internal dynamics. 

Key words. Galaxies: clusters: individual: Abell 959 - Galaxies: clusters: general - Galaxies: distances and redshifts 

1. Introduction to z ~ 0.3, e.g., Giovannini et al. I1999I see also Giovannini & 

FerettiHini Feretti |2005]l. 

In the hierarchical scenario for large-scale structure forma- 
tion, clusters of galaxies are known to be unrelaxed structures. There is growing evidence of the connection between diffuse 
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c5 ■ Instead, as the cosmic time flows clusters are interested by per- radl ° emission and clustermerging basedon X-ray data (e.g., 
turbing merging processes that constitute an essential ingredi- Bohnnger & Schuecker [2002l Buote IMS. Studies based on 
ent of their evolution. Much progress has been made in the past a lar 8 e number of clusters have found a significant relation be- 
decade in the observations of the signatures of merging pro- /ween the radi o and the X-ray surface brightness (Govom et al. 
cesses (see Feretti et al. 1200251 for a general review). A recent l™aj IZQOlbJl and connections between the presence of radio- 
aspect of these investigations is the possible connection of clus- halos/rehcs and irregular and bimodal X-ray surface brightness 
ter mergers with the presence of extended, diffuse radio sources: distribution (Schuecker et al.| 2001| >. 

halos and relics. The synchrotron radio emission of these sources Optical data are a powerful way to investigate the presence 

demonstrates the existence of large-scale cluster magnetic fields and the dynamics of cluster mergers (e.g., Girardi & Biviano 

and of widespread relativistic particles. It has been suggested 2002), too. The spatial and kinematical analysis of member 

that cluster mergers provide the large amount of energy neces- galaxies allow us to detect and measure the amount of substruc- 

sary for electron reacceleration up to relativistic energies and tures and to identify and analyze possible pre-merging clumps 

for magnetic field amplification (Feretti 1999; Feretti 2002a; or merger remnants. This optical information is really comple- 

Sarazin r2002l l. However, the precise radio halos/relics formation mentary to X-ray information since galaxies and intracluster 

scenario is still being debated since the diffuse radio sources are medium react on different time scales during a merger (see, e.g., 

quite uncommon and only recently can we study these phenom- numerical simulations by Roettiger et al. 119971 ). In this context 

ena on the basis of sufficient statistics (few dozen clusters up we are conducting an intensive observational and data analysis 

program to study the internal dynamics of radio clusters by using 

Send offprint requests to: W. Boschin, e-mail: boschin@tng . iac . es member galaxies. Our program concerns both massive clusters, 
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where diffuse radio emissions are more frequently found (e.g., 
Girardi et al. 2008 and refs. therein), and low-mass galaxy sys- 
tems (Boschin et al. 12008?!) . 

During our observational program we conducted an intensive 
study of the cluster Abell 959 (hereafter A959). It is a fairly rich, 
X-ray luminous, hot and massive Abell cluster: Abell richness 
class = 1 (Abell et al. [19891 ; L x (0.1-2.4 keV)=13.25xl0 44 
erg s (Popesso et al. 2004, but correcting for our < z >, 
see be low); T x = 6.95+[g keV (ASCA data, Mushotzky & 
Scharf |T9971 >; a projected 2D aperture mass within a radius of 
440"M ap -1-2 xlO 15 h'l M Q (Dahle et al. 120021 in their 
cosmology). 

Optically, the cluster center is not dominated by any sin- 
gle galaxy - it is classified as Bautz-Morgan class III (Abell 
et al. 1989) - but it has a core region consisting of many early- 
type galaxies of similar brightness (Dahle et al. 2002). The gen- 
eral SE-NW elongation of the A959 cluster galaxy distribution 
is matched in X-rays (ROSAT/PSPC data; Dahle et al. 120031 . 
Perpendicularly, the NE-SW direction is defined by a clear NE 
external extension and a likely SW extension in the X-ray emis- 
sion (Dahle et al. 120031 1 and by a possible very central bimodal 
structure (Ota & Mitsuda |2"0"04l . 

A highly significant mass peak (with some evidence of sub- 
structure) is seen in the weak lensing mass map (Dahle et al. 
2002 1 2003), and the mass distribution again indicates two direc- 
tions for the cluster elongation (SE-NW and NE-SW) roughly 
resembling the galaxy light distribution, although the agreement 
is not perfect in the central region (see Fig. Q] for a multiwave- 
length view of this cluster). Moreover, there is a dark matter con- 
centration at ~6' SW of the cluster center, which does not corre- 
spond to any peak in the luminosity distribution and is connected 
by a mass bridge with A959 (WL 1017.3+5931). To date the na- 
ture of this object is unknown and WL 1017.3+5931 remains a 
good candidate for an optically "dark" subcluster (Dahle et al. 
120031 . 

The presence of a radio extended emission in A959 has been 
reported by Cooray et al. dl998l l. Owen et al. (1 19991 1 lists A959 
among clusters showing a detectable diffuse radio source with 
a flux density of 3 mJy (at 1.4 GHz) and a size of ~3'(~0.8 
/i^J Mpc in our cosmology). The measurements of this prelimi- 
nary study are rough and more observations of this radio source 
are needed to clarify the possible connection with the cluster. 
However, we note that this source shows the typical sizes of 
radio halos and radio relics/gischt (see Kempner et al. 2004; 
Ferrari et al. 120081) . Moreover, its radio power (~lxl0 24 W/Hz) 
closely matches the relation logP(1.4 GHz) - logTx found by 
Cassano, Brunetti & Setti ( 120061 see their Fig. 3b) for radio ha- 
los. 

Poor spectroscopical data exist for A959. The value usually 
quoted in the literature for the cluster redshift is z = 0.353, 
but Irgens et al. (12002ft determined z = 0.286 using four galax- 
ies. We carried out spectroscopic observations at the Telescopio 
Nazionale Galileo (TNG) giving new redshift data for 107 galax- 
ies in the field of A959 and finding z = 0.288. 

This paper is organized as follows. We present our new spec- 
troscopical data and the cluster catalog in Sect. 2. We present our 
results about the cluster structure in Sect. 3. We briefly discuss 
our results and give our conclusions in Sect. 4. 

Unless otherwise stated, we give errors at the 68% con- 
fidence level (hereafter c.l.). Throughout this paper, we use 
Hq = 70 km s _1 Mpc 1 in a flat cosmology with Qq = 0.3 and 



1 please visit the web site of the DARC (Dynamical Analysis of 
Radio Clusters) project: http://adlibitum.oat.ts.astro.it/girardi/darc 



Q.A = 0.7. In the adopted cosmology, 1' corresponds to ~ 260 
h^Q kpc at the cluster redshift. 

2. New data and galaxy catalog 

Multi-object spectroscopic observations of A959 were car- 
ried out at the TNG telescope in December 2006. We used 
DOLOPvES/MOS with the LR-B Grism 1, yielding a disper- 
sion of 187 A/mm, and the Loral CCD of 2048 x 2048 pixels 
(pixel size of 15 fan). This combination of grating and detec- 
tor results in a dispersion of 2.8 A/pix. We observed four MOS 
masks for a total of 130 slits. We acquired three exposures of 
1800 s for each mask. Wavelength calibration was performed 
using helium-argon lamps. Reduction of spectroscopic data was 
carried out with the IRAF0 package. 

Radial velocities were determined using the cross- 
correlation technique (Tonry & Davis 1979) implemented in the 
RVSAO package (developed at the Smithsonian Astrophysical 
Observatory Telescope Data Center). Each spectrum was corre- 
lated against six templates for a variety of galaxy spectral types: 
E, SO, Sa, Sb, Sc, Ir (Kennicutt [19921 . The template producing 
the highest value of i.e., the parameter given by RVSAO and 
related to the signal-to-noise of the correlation peak, was cho- 
sen. Moreover, all spectra and their best correlation functions 
were examined visually to verify the redshift determination. The 
median value of % of our successfully measured galaxy redshifts 
is ~ 7.5. In 4 cases (IDs 76, 79, 84, and 90; see TableQ} we took 
the EMSAO redshift as a reliable estimate of the redshift. Our 
spectroscopic survey in the field of A959 consists of 1 19 spectra 
with a median nominal error on cz of 56 km s _1 . The nominal 
errors as given by the cross-correlation are known to be smaller 
than the true errors (e.g., Malumuth et al. 1992; Bardelli et al. 
[T9941 Ellingson & Yee fT994l Quintana et al. |2000l . Double red- 
shift determinations for the same galaxy allowed us to estimate 
real intrinsic errors in data of the same quality taken with the 
same instrument (Barrena et al. |2007al |2007bl > . Here we applied 
a similar correction to our nominal errors; i.e., hereafter we as- 
sume that true errors are larger than nominal cross-correlation 
errors by a factor 1 .4. Thus the median error on cz is ~ 78 km s . 

We also used public photometric data from the Sloan Digital 
Sky Survey (SDSS, Data Release 6). In particular, we used r', 
i', z' magnitudes, already corrected for the Galactic extinction, 
and considered galaxies within a radius of 30' from the cluster 
center. Our spectroscopic sample is ~70% complete down to r' - 
20 within 4' from the cluster center. 

Table[T]lists the velocity catalog (see also Fig. |2]i: identifica- 
tion number of each galaxy, ID (Col. 1); right ascension and dec- 
lination, a and 8 (J2000, Col. 2); SDSS r' magnitudes (Col. 3); 
heliocentric radial velocities, v = czq (Col. 4) with errors, Av 
(Col. 5). 

3. Analysis and results 

3. 1 . Member selection and global properties 

To select cluster members out of 107 galaxies having redshifts, 
we follow a two-steps procedure. First, we perform the ID 
adaptive-kernel method (hereafter DEDICA, Pisani 119931 and 
[19961 see also Fadda et al. [T996l Girardi et al. [T996l Girardi & 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Multiwavelength image of A959 (north at the top and east to the left). The contour levels (thin levels) of a smoothed archival 
ROSAT/PSPC X-ray image are superimposed on the SDSS r' image. The peaks (thick contour levels) of the projected mass density 
field (adapted from Dahle et al. I20031 > reconstructed from the weak lensing analysis are shown. Contours recovered from the galaxy 
distributions of the three groups VI, V2, and V3 we detect analyzing the velocity field (see text) are also shown (dashed thick 
contours - one contour for each group). The two black circles highlight the first and the second dominant galaxies (I and II, see 
text). Labels indicate individual subclumps discussed in the text. 



Mezzetti 2001). We search for significant peaks in the velocity 
distribution at >99% c.l.. This procedure detects A959 as a peak 
at z ~ 0.288 populated by 94 galaxies considered as candidate 
cluster members (see Fig. |3). Out of 13 non members, nine and 
four are foreground and background galaxies, respectively. 

All the galaxies assigned to the A959 peak are analyzed in 
the second step, that uses the combination of position and ve- 
locity information: the "shifting gapper" method by Fadda et al. 
( II 9961 1 . This procedure rejects galaxies that are too far in veloc- 
ity from the main body of galaxies and within a fixed bin that 
shifts along the distance from the cluster center. The procedure 
is iterated until the number of cluster members converges to a 
stable value. Following Fadda et al. ( 1 19961 1 we use a gap of 1000 
km s _1 - in the cluster rest frame - and a bin of 0.6 Mpc, or 
large enough to include 15 galaxies. 



The choice of the cluster center is not obvious. No evident 
dominant galaxy is present, but rather there are two luminous 
galaxies (IDs. 53 and 82, hereafter I and II dominant galax- 
ies) defining the SE-NW direction already shown by X-rays 
(Dahle et al. 120031 1. Moreover, the gravitational lensing analysis 
brings out several mass peaks (Dahle et al. 2003). As the clus- 
ter center, hereafter we assume the position of the peak of X-ray 
emission as listed by Popesso et al. d2004t [R.A.=10 h 17 m 35 ?04, 
Dec.=+59°33'27 '!! (J2000.0)] which is very close to the I dom- 
inant galaxy. After the "shifting gapper" procedure, we obtain a 
sample of 81 fiducial cluster members (see Fig.|4|i. 

The 2D galaxy distribution analyzed through the 2D 
DEDICA method only shows one peak [at R.A.=10 h 17 m 36 ?3, 
Dec.=+59°34'09 '.'0 (J2000.0)], which is very close to the I dom- 
inant galaxy (ID. 53), too. Using this alternative cluster center, 
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Fig. 2. SDSS r'-band image of A959 (north at the top and east to the left). Circles and boxes indicate cluster members and non- 
member galaxies, respectively (see Table[TJ. 



we verify the robustness of our member selection obtaining the 
same sample of fiducial members as above. 

By applying the biweight estimator to the 8 1 cluster mem- 
bers (Beers et al. 119901 1, we compute a mean cluster redshift of 
(z) = 0.2883+ 0.0004, i.e. (v) = 86442+130 km s _1 . We es- 
timate the LOS velocity dispersion, cry, by using the biweight 
estimator and applying the cosmological correction and the stan- 
dard correction for velocity errors (Danese et al. 1980). We ob- 
tain cr v = 1 170^3 km s -1 , where errors are estimated through a 
bootstrap technique. 

To evaluate the robustness of the cr v estimate, we analyze 
the velocity dispersion profile (Fig. [5}. The integral profile is 
roughly flat suggesting that the value of cr v is quite robust. 
However we notice a slight increasing in the external cluster 
regions confirmed by the behavior of the differential velocity 
dispersion profile that increases, rather than decreases, going to 
larger and larger radii. The <x v enhancement is likely due to the 
NE structure we detect and discuss in the following section. A 
more reliable estimate of the global cry then seems the value 
within 0.8 Mpc where cr v = 1025+^ 4 km s _1 . 

Making the usual assumptions (cluster sphericity, dynami- 
cal equilibrium and galaxy distribution tracing the mass distri- 



bution), we can compute global virial quantities using our esti- 
mate of <t v . Following the prescriptions of Girardi & Mezzetti 
(2001), we assume for the radius of the quasi-virialized region 
R VlI = 0.17 x cr v /H(z) - 2.15 Mpc - see their Eq. 1 after in- 
troducing the scaling with H(z) (see also Eq. 8 of Carlberg et al. 
|19971 for/?7nn). Thus the cluster is sampled out to a large fraction 
of Ruoq. We compute the virial mass (Limber & Mathews 1960; 
see also, e.g., Girardi et al. |1998l >: 

M = 37r/2 ■ (TyRpv/G - SPT, (1) 

where SPT is the surface pressure term correction (The & White 
1986), and Rpy is a projected radius (equal to two times the pro- 
jected harmonic radius). The value of Rpy depends on the size of 
the sampled region and possibly on the quality of the spatial sam- 
pling (e.g., whether the cluster is uniformly sampled or not). On 
the whole cluster region sampled out to 1 .48 hzl Mpc, we com- 
pute R PV (R < 1 .48 h'l Mpc) = ( 1 .24 + 0.09) Mpc. The value 
of SPT strongly depends on the radial component of the veloc- 
ity dispersion at the radius of the sampled region and could be 
obtained by analyzing the velocity-dispersion profile, although 
this procedure would require several hundred galaxies. By com- 
bining data on many clusters sampled out to about R vll a typical 
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Table 1. Velocity catalog of 107 spectroscopically measured 
galaxies in the field of A959. 
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IDs in italics indicate non-cluster galaxies. Galaxy IDs. 53 and 82 
(in boldface) indicate the first and the second dominant galaxies of 
the cluster (I and II; see text), respectively. 



Table 1. Continued. 
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b A visual inspection of the SDSS r' image suggests that the r' mag- 
nitude of galaxy ID 69 listed in the SDSS catalog is incorrect, possi- 
bly due to the contamination of a very close galaxy. We recompute 
the r' of ID 69 and present our estimate here. 



val ue of SPT~ 20% is obtained (Carlberg et al. [T997t Girardi et 
al. [1998). We assume the same correction for A959. Thus, we 
obtain M(< R = 1.48 fc^Mpc) = 1.15^ xlO 15 h~l M Q . 

Since A959 is not sampled out to R v { r and the sample is 
not magnitude complete we decide to use an alternative esti- 
mate of Rpy based on the knowledge of the galaxy distribution. 
Following Girardi et al. 1998 (see also Girardi & Mezzetti lzTjOl] ), 
we assume a King-like distribution with parameters typical of 
nearby/medium-redshift clusters: a core radius R c = 1 /20 x R w \ r 
and a slope-parameter /3g t = 0.8; i.e., the galaxy volume den- 
sity at large radii goes as r~ 3/3/ " = r~ 2A . We obtain Rpy(R < 
2.15 /z^jMpc) = 1.60 hz\ Mpc, where a 25% error is expected 
(Girardi et al. [1998b . We obtain M(< R wil = 2.15 h^Mpc) = 
1-47^43 xl ° 15 h ia M e' which rescaled using the above model 
gives M(< R = 1.48 h^Mpc) = 1.10 xlO 15 h~\ M Q , in good 
agreement with the direct mass estimate given in the previous 
paragraph. 



3.2. Velocity distribution 

We analyze the velocity distribution to look for possible devia- 
tions from Gaussianity that might provide important signatures 
of complex dynamics. For the following tests the null hypothesis 
is that the velocity distribution is a single Gaussian. 

We estimate three shape estimators, i.e. the kurtosis, the 
skewness, and the scaled tail index (see, e.g., Beers et al. 119911 1. 
According to the value of the kurtosis (2.319), the velocity dis- 
tribution is light-tailed and marginally differs from a Gaussian 
at the 90 - 95% c.l. (see Table 2 of Bird & Beers[l993l). 

Then we investigate the presence of gaps in the velocity dis- 
tribution. A weighted gap in the space of the ordered velocities 
is defined as the difference between two contiguous velocities, 
weighted by the location of these velocities with respect to the 
middle of the data. We obtain values for these gaps relative to 
their average size, precisely the midmean of the weighted-gap 
distribution. We look for normalized gaps larger than 2.25, since 
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in random draws of a Gaussian distribution they arise at most in 
about 3% of the cases, independent of the sample size (Wainer 
and Schacht 119781 see also Beers et al. |1991| l. We detect two 
significant gaps (at the 97% c.l.) that separate the main cluster 
in three groups of 20, 37 and 24 galaxies (see Fig. |6), hereafter 
named VI, V2, and V3, respectively (see Table|2]for their main 
kinematical properties). 

When compared two by two through the 2D Kolmogorov- 
Smirnov test (Fasano & Franceschini 1987 ), these groups differ 
in spatial distribution: the VI group marginally differs both from 
the V2 and V3 groups (at the 90% and 94% c.l., respectively); 
the V2 group differs from the V3 group (at the 98% c.l.). In par- 
ticular, the V2 group forms a dense elongated group around the 
cluster center and hosts both dominant galaxies. The V3 galaxies 
are mostly located in the northeastern part of the sampled region. 
For each of the these groups, Fig. [7] shows the spatial distribu- 
tion of galaxies. The contours of the three galaxy distributions 
as recovered using the 2D DEDICA method (one for each dis- 
tribution) are shown in Fig. Q] The VI and V2 groups define 
the SE-NW direction and, in particular, the V2 group shows a 
bimodal 2D distribution where the two peaks correspond to the 
two dominant galaxies. The V3 group defines the NE-SW direc- 
tion. 
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Fig. 3. Redshift galaxy distribution. The solid line histogram 
refers to (94) galaxies assigned to the cluster according to the 
DEDICA reconstruction method. 

Table 2. Kinematical properties of galaxy subclumps 



Clump 


a, S (J2000) N g 


< v > 


< 




density peak center 


km s -1 


km s 1 



Vl b 

V2 b 

V3 b 

KMM1 

KMM2 

KMM3 

DS -NE 



10 17 21.58, 
10 17 35.08, 
10 17 46.07, 
10 17 40.87, 
10 17 34.46, 
10 18 03.55, 
10 18 05.90, 



+59 35 10.3 
+59 34 17.4 
+59 34 20.0 
+59 32 35.2 
+59 34 24.7 
+59 36 15.7 
+59 36 14.6 



20 84611 +98 

37 86295 + 63 

24 88229 + 103 

24 85464 + 174 

40 86351 + 153 

17 88241 ± 190 

6 88945 + 319 



423 + '° 3 
378+5 

4Q 1 +66 

831 + 'f 
954!f 

7C0 + I6O 
' JJ -108 

65 7 +727 



We use the biweight estimator and the gapper estimator by Beers et 
al. (1990) for samples with N e > or < 15, respectively (see Girardi 
etal. [T99ll >. 

The estimate of o" v should be considered a lower limit in these sam- 
ples since the membership assignment might lead to an artificial 
truncation of the tails of the velocity distributions (e.g., Bird ! 1994t . 



Finally, we notice that the luminous galaxies are preferen- 
tially located in the V2 group tracing, in particular, the cluster 
core. When comparing the magnitude distribution of V2 galax- 
ies with that of V1+V3 galaxies, we find a difference at the 98% 
c.l. according to the ID Kolmogorov-Smirnov test (see, e.g., 
Ledermann 19821. 



3.3. "3D" substructure analysis 

The existence of correlations between positions and velocities of 
cluster galaxies is a footprint of real substructures. Here we use 
different approaches to analyze the structure of A959 combining 
velocity and position information. 
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Fig. 4. Rest frame velocity vs. projected clustercentric distance 
for the 94 galaxies in the main peak (Fig. [3]l showing galax- 
ies detected as interlopers by our "shifting gapper" procedure 
(crosses). Solid squares indicate cluster members detected as the 
northeastern substructure in the Dressler-Shectman test (DS- 
NE). 



We analyze the presence of a velocity gradient performing a 
multiple linear regression fit to the observed velocities with re- 
spect to the galaxy positions in the plane of the sky (e.g, Boschin 
et al. 120041 and refs. therein). We find a position angle on the 
celestial sphere of PA = 85^ degrees (measured counter- 
clockwise from north), i.e. higher-velocity galaxies lie in the 
western region of the cluster (see Fig. [8]). To assess the signif- 
icance of this velocity gradient, we perform 1000 Monte Carlo 
simulations by randomly shuffling the galaxy velocities, and for 
each simulation we determine the coefficient of multiple deter- 
mination (RC 2 , see e.g., NAG Fortran Workstation Handbook 
1986). We define the significance of the velocity gradient as 
the fraction of times in which the RC 2 of the simulated data is 
smaller than the observed RC 2 . We find that the velocity gradient 
is significant at the 99.6% c.l.. 

We combine galaxy velocity and position information to 
compute the A-statistics devised by Dressier & Shectman 
(1988). This test is sensitive to spatially compact subsystems 
that have either an average velocity that differs from the cluster 
mean or a velocity dispersion that differs from the global one, or 
both. We find A = 126 for the value of the parameter that gives 
the cumulative deviation. This value is a significant indication of 
substructure (at the 99.9% c.l.) as assessed by computing 1000 
Monte Carlo simulations, randomly shuffling the galaxy veloc- 
ities. Figure [9] shows the distribution on the sky of all galaxies, 
each marked by a circle: the larger the circle, the larger the devia- 
tion 5j of the local parameters from the global cluster parameters, 
i.e. the greater the evidence for substructure. 

To better point out galaxies belonging to substructures, we 
resort to the technique developed by Biviano et al. (2002 see 
also Boschin et al. 2006; Girardi et al. 2006), who used the in- 
dividual values of the Dressler-Shectman method. The crit- 
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Fig. 5. Differential (big circles) and integral (small points) pro- 
files of mean velocity (upper panel) and LOS velocity disper- 
sion (lower panel). For the differential profiles we plot the val- 
ues for five annuli from the center of the cluster, each of 0.3 
Mpc (large symbols). For the integral profiles, the mean and 
dispersion at a given (projected) radius from the clump-center is 
estimated by considering all galaxies within that radius - the first 
value computed on the five galaxies closest to the center. The er- 
ror bands at the 68% c.l. are also shown. In the lower panel, the 
horizontal line represents the X-ray temperature with the respec- 
tive 90 per cent errors (Mushotzky & Scharf fT9971 l transformed 
in cry assuming the density-energy equipartition between gas 
and galaxies, i.e. j0 spe c = 1 (see text). 



ical point is to determine the value of (5, that optimally indi- 
cates galaxies belonging to substructure. To this aim we consider 
the (5,— values of all 1000 Monte Carlo simulations used above. 
The resulting distribution of 5, is compared with the observed 
one finding a difference at the 99.6% c.l. according to the ID 
Kolmogorov-Smirnov test. The "simulated" distribution is nor- 
malized to produce the observed number of galaxies and com- 
pared with the observed distribution in Fig.[T0j the latter shows a 
tail at high values. The inspection of Fig.[10]suggests that galax- 
ies with S, > 3.5 likely are in substructures. These galaxies are 
indicated with heavy circles in Fig. [9] The six galaxies located 
in the northeastern region of the cluster are likely to form a sub- 
clump, hereafter DS-NE (see Tableland Fig.@]l. 

We also use here the Kaye's mixture model (KMM) test to 
find a possible group partition of the velocity distribution (as 
implemented by Ashman et al. 1 19941 >. The KMM algorithm fits 
a user-specified number of Gaussian distributions to a dataset 
and assesses the improvement of that fit over a single Gaussian. 
In addition, it provides the maximum-likelihood estimate of the 
unknown n-mode Gaussians and an assignment of objects into 
groups. This algorithm is usually used to analyze the velocity 
distribution where theoretical and/or empirical arguments indi- 
cate that the Gaussian model is reasonable. Here the ID KMM 
test fails in confirming the three groups detected by the weighted 
gap method. However, since in this case there is a significant cor- 
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Fig. 6. Upper panel: rest frame velocity histogram for the 81 
cluster members. The bigger and smaller arrows indicate the ve- 
locities of the I and II dominant galaxies. The three Gaussians 
correspond to the three groups (VI, V2, and V3, from left to 
right) found in the velocity distribution (see Table 2). Lower 
panel: stripe density plot where the arrows indicate the positions 
of the significant gaps. 
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Fig. 7. Spatial distribution on the sky of the 8 1 member galax- 
ies showing the groups recovered by the weighted gap analysis. 
Open squares, open circles and solid circles indicate VI, V2 and 
V3 groups. The larger and smaller crosses indicate the positions 
of the I and II dominant galaxies, respectively. 



Fig. 8. Spatial distribution on the sky of 8 1 cluster members. The 
larger the symbol, the smaller the radial velocity. The plot is cen- 
tered on the cluster center. The solid and faint lines indicate the 
position angle of the cluster gradient and relative errors, respec- 
tively. 
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Fig. 9. Spatial distribution of the 81 cluster members, each 
marked by a circle: the larger the circle, the larger the deviation 
5i of the local parameters from the global cluster parameters; i.e., 
there is more evidence of substructure according to the Dressier 
& Shectman test (see text). The boldface circles indicate those 
with 5i > 3.5 (see text). Out of these, the six northern galax- 
ies form the DS-NE clump. The plot is centered on the cluster 
center. 
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Fig. 10. The distribution of <5,- deviations of the Dressler- 
Shectman analysis for the 81 member galaxies. The solid line 
represents the observations, the dashed line the distribution for 
the galaxies of simulated clusters, normalized to the observed 
number. 
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Fig. 11. Spatial distribution on the sky of the 81 galaxies of the 
whole cluster showing the groups recovered by the weighted gap 
analysis. Open squares, open circles, and solid circles indicate 
KMM1, KMM2, and KMM3 groups. The larger and smaller 
crosses indicate the positions of the I and II dominant galaxies, 
respectively. 



relation between galaxy velocities and positions we have decided 
to use KMM as a 3D diagnostic. We use the three-group galaxy 
assignment given by the weighted gap analysis to determine the 
first guess when fitting three groups. The algorithm fits a three- 
group partition at the 92% c.l. according to the likelihood ratio 
test (hereafter KMM1, KMM2, KMM3 groups of 24, 40, and 17 
galaxies from low to high mean velocities). The results for the 
three groups are shown in Table [2] and Fig. QT| Although there 
are some differences between the group memberships recovered 
with the two different methods (cf. Fig. [TT] with Fig. [7]), there 
are important common features: the NE region is populated by 
the high-velocity galaxies while the cluster core, elongated in 
the SE-NW direction, is formed by the medium velocity galax- 
ies. The main difference is that the SE-NW bimodal structure is 
now somehow split, with KMM1 and KMM2 tracing the SE and 
NW peaks better, respectively. However, note that the two dom- 
inant galaxies are both contained in the medium velocity group 
KMM2. Moreover, the galaxies in SW region are now assigned 
primarily to the low velocity group. 

3.4. 2D galaxy distribution 

When applying the DEDICA method to the 2D distribution 
of the 81 galaxy members, we find only one significant peak. 
However, our spectroscopic data suffer from magnitude incom- 
pleteness, and the sampled region does not cover the whole virial 
region. To overcome these limits we use the photometric catalog 
extracted from the SDSS (Data Release 6). 

In the SDSS photometric catalog, we select likely members 
on the base of the (r'-i') - (i'-z') color-color plane (see e.g. 
Goto et al. |2002t Boschin et al. [2008 ). Out of our photometric 
catalog, we consider galaxies (here objects with r' phototype pa- 
rameter = 3) lying within +0.08 mag from the median values of 



r'-/'=0.49 and /'-z'=0.35 colors of the spectroscopically clus- 
ter members (see Fig. [12) . The value of 0.08 mag is two times 
the typical scatter reported by Goto et al. (2002) for the corre- 
sponding color-magnitude relations r'-i' vs. r' and i'-z' vs. r' . 
By adopting this alternative selection method in a region within 
30' (almost 4 R vn ) from the cluster center we select a sample of 
likely cluster members with r' < 2 1 . The result of the application 
of DEDICA to this sample is shown in Fig.Q~3] The galaxy dis- 
tribution shows a shape elongated along the SE-NW direction in 
the central regions. The most significant peaks (at the > 99.9% 
c.l.) correspond to the I dominant galaxy, the II dominant galaxy, 
a western peak well outside of R wn (hereafter "W" peak), the NE 
peak at <5'and a far NE peak at ~18' from the cluster center 
(hereafter "far NE" peak). 

To probe the robustness of these detections, we also apply 
the Voronoi tessellation and percolation (VTP) technique (see 
Ramella et al. 1200 It Barrena et al. 2005). This technique is non- 
parametric and does not smooth the data. As a consequence, it 
identifies galaxy structures irrespective of their shapes. In partic- 
ular, we run VTP in a region within 30' from the cluster center 
selecting galaxies with r' < 22 and the same color constraints as 
above. The result of the application of VTP is shown in Fig. [14] 
VTP is run 4 times adopting four detection thresholds: galaxies 
identified as belonging to structures at the 95%, 98%, 99%, and 
99.9% c.ls. VTP reveals the elongation of A959 in the SE-NW 
direction and confirms the NE peak, the W, and far NE peaks. 

4. Discussion and conclusions 

Our analysis based on member galaxies shows that A959 is a 
strongly substructured cluster. In particular, we detect an NE 
clump, at ~5' from the cluster center, using the weighted gap 
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Fig. 12. i'—zf vs. r'-i' diagram for galaxies with available spec- 
troscopy is shown by small crosses. The faint box is centered 
on the median value for colors of member galaxies (circles) and 
encloses galaxies having colors in the range of +0.08 mag from 
median values. 
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Fig. 13. Spatial distribution on the sky and relative isodensity 
contour map of likely cluster members (according to the selec- 
tion in the color-color plane) with r' < 21, obtained with the 
DEDICA method. Dominant galaxies I and II, the NE, W, and 
far NE peaks are indicated. The circle indicates the likely virial- 
ized region. 



analysis (V3 group), the Dressler-Shectman statistics (DS-NE 
group), the 3D KMM analysis (KMM3 group) and the 2D anal- 
ysis of the galaxy distribution. This clump has a rest frame veloc- 
ity of Vrf ~ +1900 km s _1 and is the likely cause of the velocity 
gradient we observe in A959. The NE clump is also detected in 
ROSAT/PSPC X-ray data and in the lensing gravitational map 
by Dahle et al. (2003). This clump is also aligned with the SW 
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Fig. 14. Galaxies belonging to structures as detected by the 
Voronoi tessellation and percolation technique. The algorithm 
is run on a sample of galaxies with r' < 22 selected in the color- 
color plane (see text) within 30' from the center of A959. Open 
squares, crosses, asterisks, and solid circles indicate galaxies in 
structures at the 95%, 98%, 99%, and 99.9% c.ls., respectively. 
The NE clump and the W and far NE groups are indicated. 



mass concentration WL 1017.3+5931 whose nature is very in- 
triguing. 

The central cluster region is not well relaxed, too. We find ev- 
idence of a central, dense structure elongated along the SE-NW 
direction as shown by the the weighted gap analysis (V2 group), 
the 3D KMM analysis (KMM2 group) and the 2D analysis of 
the galaxy distribution. The SE-NW central structure has an in- 
termediate velocity and hosts more luminous cluster galaxies; 
thus it likely represents the main part of the cluster still showing 
the trace of a past merger. The apparent elongation of this struc- 
ture is connected with the presence of the two dominant galaxies 
and some galaxies surrounding them, as shown by the two peaks 
detected by DEDICA and VTP in the 2D galaxy distributions. 
Close couples of dominant galaxies having different velocities 
are often observed in clusters (Boschin et al. 120061 Barrena et 
al. 12007 a) and are the likely tracers of a previous cluster merger. 
Indeed cluster merger is thought to be the cause of the formation 
of dumbbell galaxies (e.g., Beers et al. 1992; Flores et al. 2000). 
With respect to other cases we have studied (e.g., clusters Abell 
773 and Abell 2744; see Barrena et al. |2007al and Boschin et 
al. 1200 6, respectively), the two dominant galaxies are closer in 
the velocity space and more distant in sky positions suggesting 
that the direction of the old, past merger is almost parallel to the 
plane of the sky. 

The SE-NW elongation is also shown by gravitational lens- 
ing and ROSAT X-ray data (Dahle et al. l200lT l. The ellipticity of 
the X-ray surface brightness appears less pronounced than those 
of mass and galaxy distributions, in agreement with the colli- 
sional component of the cluster having already had some time to 
relax. 
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Table 3. Properties of structures not listed in Table [2] 
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Structure 


Q-,(5(J2000) 


Source 


Additional information 



WL 1017.3 + 5931 
NE extension 
WSW 
W 

far NE 



10 17.3, +59 31 gravitational lensing (Dahle et al. l2003l 

10 18.2, +59 37 X-ray (ROSAT archival data) 

10 16.7, +59 32 X-ray (ROSAT archival data) 

10 15.9, +59 35 2D gal. distr. (SDSS, Koester et al. 2007) 

10 19.2, +59 46 2D gal. distr. (SDSS, Koester et al. 2007) 



no optical counterpart 

likely connected with our high velocity clump V3 
likely infalling group 

MaxBCG J153.93477 + 59.57870 (z phot ~ 0.281) 
MaxBCG J 154.78308 + 59.76784 (z phot ~ 0.284) 



The above discussion suggests that A959 is forming along 
two main directions of mass accretion (SE-NW and NE-SW). 
Indeed, the region around this cluster is quite complex and de- 
serves a brief discussion. By running the VTP technique and 
DEDICA on the 2D galaxy distribution we find two significant 
galaxy concentrations: a group at ~13' W of A959 and a far 
NE one at ~18' from the core of the cluster. These two groups 
were first detected by Koester et al. d2007) using the maxBCG 
red-sequence method in SDSS photometric data. They com- 
puted photometric redshifts of ~0.28 1 and ~0.284 for the W and 
NE groups, respectively. Thus, they are at the same redshift of 
A959, at projected distances from it of -3.4 Mpc (W group) 
and ~4.7 Mpc (NE group), possibly representing the most 
important knots of the filamentary large-scale structure around 
the cluster, after A959 itself. Closer to A959 (at ~6' WSW 
from the cluster center), ROSAT archival data reveal a diffuse 
X-ray emission probably associated to a small galaxy group. 
Unfortunately, we have no redshift in this region, but at the cen- 
ter of the X-ray emission there is a brilliant galaxy (only slightly 
less luminous than the I and II dominant galaxies) showing the 
typical colors of the early-type galaxies of A959. This galaxy 
(known as NVSS J101642+59222 2) is a n evident pointlike ra- 
dio source in NVSS (Condon et al. 1 19981 1 archival data, too. We 
suspect that this radio galaxy is at the center of a galaxy group 
falling in the gravitational well of A959. In Table [3] we summa- 
rize the properties of the WSW group, the W group and other 
structures not listed in Table [2] 

As for the global properties of the cluster, the value we find 
for the velocity dispersion of all A959 cluster members is very 
high, <t v = 1170^73 km s _1 . When considering a more central 
region to avoid the NE subclump we still obtain a high value 
cr v (< R = 0.8 /i^Mpc) = 1025!'° 4 km s" 1 . This value is in good 
agreement with the average X-ray temperature (Mushotzky & 
Scharf 1 19971 1 when assuming the energy density equipartition 
between gas and galaxies (i.e. yS spe c = 10, see Figj5j and leads 
to a virial mass of M(< R = 1.48 /z^'Mpc) = 1 . 15^" j| 
xlO 15 M Q . We compare our mass estimate with that de- 
rived from the gravitational lensing analysis. Dahle et al. (2002) 
obtained a projected mass M pro j(< 1.9 /i^Mpc) = 1.6 - 3.2 
xlO 15 h'l M (in our cosmology). To rescale and project our 
mass, we assume two alternative mass distributions: a King- 
like mass distribution (see above) or an NFW profile where 
the mass-dependent concentration parameter c is taken from 
Navarro et al. (1997) and rescaled by the factor 1 + z (Bullock 
et al. 1200 It Dolag et al. 12004ft . Our projected mass estimate is 
M proj (< 1.9 /z^Mpc) = (1.5 - 1.9) xlO 15 M B , where the 
projection is done assuming that the cluster mass distribution is 
truncated at one or at two virial radii. The lensing mass is con- 



3 Apec = °"\ /(kT/fim p ) with ji = 0.58 the mean molecular weight and 
m p the proton mass. 



sistent with our projected estimate within their errors, even if 
the former is likely to be a bit higher, possibly due to the young 
dynamical state of the cluster. 

In conclusion, A959 shows the typical features of most clus- 
ters hosting halo/relic radio sources, i.e. a high mass (higher than 
0.7 x lO 15 h~l M within 2h^ Mpc; see Giovannini & Feretti 
2002) and a young dynamical state (diffuse radio sources ap- 
pear in clusters possessing substructures, i.e. experiencing large 
departures from a virialized state; see, e.g., Buote|2001). This 
makes of this cluster a good candidate for deep radio studies. 
More observations of the diffuse radio source will allow to clar- 
ify its connection with the internal dynamics of this cluster. 
Moreover, deep X-ray observations will be useful for determin- 
ing the merging phase of the NE clump. 
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